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Abstract. We develop a consistent formalism in order to explore the effects of 
density and spin fluctuations on the quasi-particle properties and on the pairing critical 
temperature of a trapped Fermi gas on the attractive side of a Feshbach resonance. We 
first analyze the quasi-particle properties of a gas due to interactions far from resonance 
(effective mass and lifetime, quasi-particle strength and effective interaction) for the 
two cases of a spherically symmetric harmonic trap and of a spherically symmetric 
infinite potential well. We then explore the effect of each of these quantities on Tc and 
point out the important role played by the discrete level structure. 
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1. Introduction 

In the last few years the progress in cooling and manipulation of atomic Fermi gases has 
been impressive PQ 12 13 IH E] ■ The most notable feature of atomic gases is the presence of 
scattering resonances (Feshbach resonances), which can be addressed using an external 
magnetic field. The resonance is induced in the scattering between two atoms in different 
internal states, typically hyperfine states, and results in the divergence of the two-body s- 
wave scattering length ap- Feshbach resonances allow the experimental study of a Fermi 
gas at various interaction regimes. The crucial feature is that, by varying the value of 
the scattering length, one can explore different kinds of fermionic superfuidity, ranging 
from the weak-coupling BCS superfluid, for ap small and negative, to the strong-couphng 
unitarity regime, ap oo, and finally to the Bose-Einstein condensation regime, ap > 0, 
where the fermions couple pair- wise to form bosonic molecular bound states. In fact, 
it has been pointed out some years ago that BCS and BEC are two aspects of one and 
the same phenomenon, and the basic theory which connects the two was then developed 
013 El- Recently, this theory came again to the attention of the scientific community 
in the context of cold atomic Fermi gases [OlCniE]- One can show that BCS and BEC 
appear at the opposite sides of the resonance if one keeps track of the mean-field normal 
and anomalous pair correlation functions. This essentially corresponds to extending 
mean-field BCS theory to the whole crossover to BEC. Although the resulting theory is 
remarkably elegant, it is only qualitatively correct. The absence of important ingredients 
from the theory as the fluctuating parts of the pair correlation functions and the higher 
order correlations leads to various flaws in the quantitative predictions. Keeping the 
full correlations is of course impracticable unless one resorts to quantum Monte Carlo 
calculations [121 CHI IHj • To pinpoint and understand the important physics contained 
in the higher order correlations, however, it is also useful to apply approximate methods. 
This is the line we follow in the present work. 

We are interested in improving our understanding of the properties of a Fermi gas 
as the interaction strength increases originating from the weak-coupling BCS regime. 
In particular the key quantity we calculate is the critical temperature for pairing Tc. 
Moving towards the resonance, higher order correlations and fluctuations over the mean- 
field become important. One may expect the emission and reabsorption of density and 
spin fluctuations to play an important role. There are very close analogies between a 
gas in these conditions and a Fermi liquid. One can build a theory along similar lines 
introducing the concept of quasi-particles displaying an effective mass, a finite lifetime 
and an effective interaction due to coupling with fluctuations. The current experiments 
with Fermi gases have one more particularly interesting aspect in that the atoms are 
usually confined in harmonic potentials. The finite size plays an important role when 
the number of particles is not large, as is also the case in optical lattices with few atoms 
per lattice site, when tunneling between neighboring sites is suppressed and atoms are 
subjected to an effective harmonic confinement. Both these situations are of interest 
and have close relations with the physics of finite nuclei and atomic clusters (cf. e. g. 
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IT^ [T7] ) . In this article we first develop a formalism appropriate to finite size systems 
confined by a general, isotropic potential. Subsequently we apply the formalism for the 
two specific cases of an infinite square well, which we compare with an homogeneous 
system, and of a spherically symmetric harmonic confinement (cf. also |18|). The 
results illustrate the central role the discrete level structure plays in determining the 
properties of these systems. Finally we relate our approach to other calculations for 
(or equivalently for the T = pairing gap) made in different fields. 

The paper in organized as follows. In Sec. |2lwe present the theoretical framework. 
We show the approximations used which lead to the emergence of the renormalized quasi- 
particle properties, namely the normal self-energy (Sec. I2.1|l . the anomalous self-energy 
and the induced interaction (Sec. 12. 2|) . and the one pole approximation (Sec. 12. 3|) . 
Finally we derive the equations used to calculate Tc (Sec. I2.4|l . including the multipolar 
expansion appropriate for an isotropic system (Sec. I2.5|l . In Sec. El which is divided 
into two main subsections, we report the results of the numerical calculations. Sec. 13.11 
contains the results for the infinite square well. We discuss the quasi-particle properties 
and compare the results obtained with those of an infinite system. We also present the 
calculated values for Tc. In Sec. 13.21 the same subjects and corresponding results for 
a harmonically trapped system are presented. Section 0] contains the conclusions. In 
[Appendix A. 1| details concerning the correlation functions and the linear response of the 
system for both the uniform and harmonically trapped case are collected. 

Before proceeding a remark on notation. In a uniform system the many-body 
interaction strength is parametrized by A = gN{0), with g being the coupling constant 
(related to the scattering length ap) and iV(0) = mkp/27i^h'^ the density of states at 
the Fermi surface for a single spin orientation. For such system the Fermi momentum 
kp is related to the uniform density n of the gas by (37r^)^/^n^/^. A similar definition 
can be introduced for a trapped gas using in place of n the density at the center of the 
interacting cloud. 

2. Theoretical tools 

We consider a Fermi gas of atoms in two internal states with equal numbers {N/2) 
described by the Hamiltonian 



Here Hq = —h^V'^/2ma + Vext{r), with nia being the atomic mass, and | and | 
conventionally label the two atomic internal states. The only assumptions made at 
this level on the external potential are that it is isotropic and that it acts in the same 
way on both atomic states. These assumptions imply that the following relation holds 
for the chemical potentials //^ = /i| = fi. Finally, we are here interested in the BCS side 
of the resonance where g < 0. 





(1) 
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The properties of the system at any given temperature T can be described 
introducing the normal and anomalous propagators ^^^(r, r', tu^) and ^^i{r,r',uJn)- 
They satisfy the Dyson equations (cf. e.g. ref. pi 



%(r,r',a;„) =^Yr,r',uj^) (2) 
+ 



j d^s d^s' "S^^ (r, s, ujr^) (s, s', (s', r', u)^ 



- Vr|j,(s, s',u;„)^/^(s',r',a;„)] 
^/^(r,r',a;„) = ^ rf^s dV ^° (s, r, -^„)[iy/^(s, s', ^„)%(s', s, (3) 

+ Su(s',s, -^„)^|^(s,r',^„)]. 

The fermionic Matsubara frequency is as usual defined by ujn = {2n + l)T[/h(3, where 
(3~^ = ksT, with T being the temperature of the system and ks the Boltzmann constant. 
Eo-o- and W^-a are the normal and anomalous self-energies respectively. Analogous 
equations to and (jlj) are satisfied by and ^-fj. The assumption of non- vanishing 
pair correlations below implies a specific choice of gauge in order to break the 
symmetry of the Hamiltonian. This is implicit in our calculations. We do not need 
to explicitly worry about the issue because we are only interested in finding the critical 
temperature Tc itself. 

In order to further proceed we now need to assume specific approximations for the 
self-energies. 

2.1. The normal self-energy 

In this work we account for the possibility for particles to emit and reabsorb density 
and spin fluctuations and adopt the following set of approximations. We assume 
E||(ci;„) = + S|^(a;„), and similarly for The first term corresponds to the 

(cu-independent) Hartree self-energy |2Uj : 

Sff(r,r') = -|^e-"^u(r,r',a;„)5(r-r') (4) 

n 

In practice it is generally not needed to consider the full in Eq. Q. It is instead 
sufficient to first solve the self-consistent Hartree problem in which is replaced by 
(with [^lY]^^ = ~ ^u)' ^^"^ then use this set of functions to calculate 

the effect of phonons. The self-consistent Hartree problem is solved by expressing the 
Hartree Green's functions as 

g-'/(r.r'..„)^^ ^:f' fef' . (5) 

The quantity (f)v is the wave function of the level v in the Hartree approximation, and 
is equal to e'l,^ — fi, with e[j:-^ being the energy of the level. In Eq. (jSj we omitted the spin 
indices since the Hartree Green's function is the same for the two species. The Hartree 
self-energy is on the other hand given by ^^(r, r') = 6{r — r') J2„ l0i^(i")P'^F(^!^)5 with 
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Figure 1. Phonon-induced self-energy (r, r', Wn) (Eqs. © and Q). The 
continuous lines correspond to Hartree Green's functions $^''^(r, r', Wn), the dashed 
lines to the bare interaction g(5(r— r'). The self-energy S^^^ is given by the contributions 

in (a). , which arises from coupling to spin phonons, is given by the diagrams in 
(b). The diagrams with an even number of bubbles are not included in (a) because the 
bare interaction acts only between atoms in different internal states. 



— 1/(6^^ + 1) being the Fermi distribution function. The chemical potential is 
obtained imposing 

N = 2J2 [ rfV|0,(r)|V(e.). (6) 

We deal with the self-energy Yf^ due to coupling with collective modes in the context 
of the Random Phase Approximation 121] . We thus suppose Yf^^ to be given by the sum 
of the Feynman graphs shown in Fig. where to each fermion line corresponds a 
and to each dashed line an interaction g5{T — T'). This part of self-energy can be written 
as the sum of two contributions t!^! and t!"^^ (cf. Fig. [T] (a) and (6) respectively), t!'^^ 
is given by 

1 9^ h 

Sj^^^(r, r', ^^n) - ^ S^tY^^' ^")nu(r, r', u;„ + Um) 

where 11^ is the ring approximation to the correlation function (for all the correlation 
functions, and their relations, used in this Section we refer the reader to Appendix A.l ). 
This can be written as the sum of a density and a spin part: 

(r, v',cOm) = \{ n;(r, r',u;J + Ul^{r,r',u;m)}, (7) 

with = 2m7T/hf3 being a bosonic Matsubara frequency. The second contribution is 
instead given by: 

2 q"^ h 
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|n;^ (r, r', + ujm) + (r, r', Un + uj.^) } . (8) 
Combining the two contributions one finds 

X |n^(r, r', Un + UJm) + SIll^ (r, r', Un + uJm) - 2nJ(r, r', + uJm) } 

-U\^{r,r',Un + uJm) - Xo(r,r',t^n + c^m)}- (9) 

In Eqs. (0) and © we have used IIo-^ = 11^^ = U^^ and we accounted for the fact 
that the single bubble diagram belongs both to S^^^ and to and is therefore doubly 
counted in a simple summation. For this reason we subtracted the quantity 211° in Eq. 
0. 

It is useful to introduce the spectral representations for the correlation functions: 

poo 7 / r) / 

U;{r,r',uJm)= — Im[n^(r,r',a;')]^— ^, (10) 

Jo + ^ 

poo J I o / 

n;^(r,r',^^)= / — Im[n;^(r,r',u;')]^— ^, (11) 

n%v,r',cOm)= —Mn%v,v',co')]^^^, (12) 
Jo n ^ ^,n + ^ 

where njj^(r, r', tu), njj.^(r, r', cj) and n°(r,r',co') are the retarded correlation functions 
(obtained using the analytic continuation cOn co + iri), and we have used Iii{r, r', a;) = 
-nKr,r',-a;). 

Inserting these relations, together with Eq. 0, in Eq. and carrying out the 
frequency summation one finds 



PC 



TT 



Im[2n'|:^(r,r',u;') - ff!f^(r, r', u;') - Xo{r,r' ,uj')] 

where UBiyj) = l/(e^^ — 1) is the Bose distribution function. 

The self-energy of the quasi-particle in the level u is given by Il^^{uJn) = 
J d'^r d^r' (f)l{T)(l)^(r')T.P'^(r,r',uJn) and therefore 



/•OO 7 / 
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Figure 2. Plionon-induced interaction l/|""'(r, r', a;„) (see Eq. H16(l ^. The diagrams 
witli an odd numbers of bubbles are not included in (a) because we are considering 
only the interaction between quasi-particles with opposite spin states. 



where 



X Im[2m. (r, r', u) - Hf , (r, r', u) - xo(r, r', u)] 



(15) 



2.2. The anomalous self- energy 

The anomalous self-energy in Eqs. (jH)) and is given by (cf. [T9] ) 

where V.^^'^ is an effective interaction between quasi-particles in states t and |. 
Accounting also for the possibility for particles to interact exchanging density and 
spin fluctuations, V.^f-^{r,r',u!m) is the sum of two terms: the (cu-independent) bare 



interaction gS{r — r') and the induced interaction 
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(16) 



which corresponds to the diagrams shown in Figure El 

ITp is related to the exchange of density modes, 11^. to the exchange of spin modes. 
Introducing Eqs. ([111)) and (fTT|) in Eq. (fT^ yields 

^2 roo^^> 2tu' 



F/f'^(r,r',^^) = ^ 



/2 



m 

xlm[n;(r,r',c.')-3n;^(r,r',a;')]. (17) 

Since the imaginary parts of the retarded correlation functions are always negative, 
we see that the density modes induce an effective attractive interaction, while the 
spin modes induce a repulsive one. In the Hartree-Fock basis the anomalous self- 
energy, as well as depends in principle on two indices ui and z/2 as Wy^^^iun) = 
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J d^rd^r' (j)i,-^{jc)(j)y^{Y')W^i{r,Y',uJn). In the same basis, Eq. (fT^ becomes: 



where 



Solving the Dyson equations with the full v dependence is rather complicated. A 
consistent simplification may be attained by keeping only the matrix elements between 
states connected by the time reversal {(py^ = and (p^^ = (pl^)- This is in keeping 
with the fact that those matrix elements are the largest because the overlap of the 
wavefunctions in Eq. (|T8|l is maximal (see for instance Ref. |22t l23]). With this 
assumption both the anomalous self-energy and the anomalous Green's function depend 
on a single index and Eq. (|T8|l becomes: 

= -IY1 Kl'M^u'{ujn + cu^), (18) 

with 

1 ref f / \ f duj' . 2cj' . , 

g^yi being the matrix element of the bare interaction and 

^ j rfWr'0,(r)0:(r')0:,(r)0,,(r')Im[n;(r,r',a;) -3n;^(r,r',c.)]. 
The Dyson equations Eq. and on the other hand take the simple form: 

%{uJn) = G^icOn) - G,{uOn)W,{u:n)^lM, (20) 

and 

= Gy{-UJr^)Wl{uJnW,{uJn), (21) 

where 

GM.)- ^ , ■ (22) 

Close to Tc one can keep only the linear contribution of the anomalous self-energy 
and obtain for the anomalous propagator 

^,{uJn)W,{yjn)G,{uJn)G,{-uJn). (23) 



2.3. One pole approximation and quasi-particle properties 

We now proceed by expanding the propagators close to the poles which correspond to 
the quasi-particle excitation energies. As will be shown below, this approximation holds 
in the weak and likely also in the medium coupling regimes. By weak coupling we mean 
A < 0.2 and by medium coupling 0.2<A<0.5. In order to carry out the one pole 
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approximation on the finite temperature propagators (the T = case is treated e. g. in 
Ref. one needs to introduce the spectral representation 

GuM = / ^ -. 7 24 

where the spectral density p is given by 

p.M-z[G:^M-G^M] (25) 

with G^[yj) = Gu{uJn)\Lj„^uj-ir, and G^{uj) = G^{ujn)\cvn^uj+iri being the advanced and 
retarded correlation functions respectively. The corresponding self-energies satisfy the 
relations 

ImEf '^(cu) = -ImSf '^(cj) > Vcj (26) 

and 

ReEf'^(u;) =ReEf'^(cj) Vcu. (27) 

If |ImE^''(u;)| <^ |^i, + ReS^^(co')|, the inverse advanced Green's function is approximately 
given by 

[Gtico)]-' ^ Z;\co - e.) - zlmSf '^(e.), (28) 

where satisfies the relation = + ReS^''(e,y). The analogous expression holds 
for the retarded Green's function. With the definition j^, = ImE^'^''^(e;^) the spectral 
function becomes 

p.{u) = 2Z,- (29) 

with the characteristic lorentzian shape of width 7,^ centered about the renormalized 
quasi-particle energy ey. The factor Zy is defined by 



1 



d{tko) 



< 1, (30) 



and it physically represents the strength of the lorentzian quasi-particle peak, as it is 
shown by the relation: 

J^P^i^')- (31) 



(32) 



Substitution into Eq. (j24|) finally gives 



GyiuJr, 



9{uJn) _^ 0{-UJ„] 



ihujn - ey + i^y ifkOn-tu-iiu_ 
(see Ref. 123)- The possibility of performing the one pole approximation relies on 
the assumption that the imaginary part of the phonon-induced self-energy is small as 
compared to the quasi-particle energy in the proximity of the Fermi surface. This is 
equivalent to requiring the quasi-particle lifetime to be long. The discussion on the 
validity of this assumption will be presented in Sect. Ill, where we will show the results 
for the self-energy calculated using Eq. ()14|). These are also important in order to 
determine in which range of the interaction parameter the description of Fermi gases in 
terms of quasi-particles is well suited. 
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Using Eq. ()32j) the anomalous propagator in Eq. (j23jl becomes: 

1 1 



2ey \-i\fvjJn 
It allows a spectral representation itself: 



with 



271 iuJr, — UJ'' 



(33) 



2ty 



27. 



27. 



7^ (a; + e,)2 + ^2j 



Using the one pole approximation for Eq. ()34|) and the symmetry of the gap function 
Wy[uj) = W^{-uj) Eq. dnni) becomes 

1 1 



2e, 



.i\hwn\ - eu + i'Ju i\hiOn\ + ty + i'^y- 
where A^, = ZyWy{ey). The gap equation now follows by inserting Eqs. (fT^ and (jHlj) 
into Eq. (HH) 

f°° duj ^ ^ 2uj 



Zyl I^yl 



2e,, 



TT 



X 



Z I hjjJn' I — Cj,/ + 27;^/ i I hjJn' | + Cjy' + ^7i. 

with ijjn' = ^^n + ^^m- Neglecting the quasi-particle width and carrying out the Matsubara 
sum one gets 



WJuJr, 



2e,/ 



X i guu' + / Vyy'{uj) 







TT 



1 



After the analytic continuation iujn — > and multiplication by Zy one finds: 



2e,/ 

X + / —Vyyi[uj) 



+ 



.C/y + |e!/'| + -^y + \ey'\ + UJ _ 
with = ZyQyyiZyi Qnd slmlkrly for Vyyi. The terms in the curly brackets constitute 
an effective interaction matrix element between quasi-particles in states v and v' . The 
second term in particular has the form of an induced interaction obtained from second 
order perturbation theory with a fermion-boson Hamiltonian of the form: 



- ft 



Vyy,{uj) dl„bl^dy>^> + C.C. 



(34) 
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^yv„lyl(uJ) being the effective particle- vibration coupling f^. In the equation a =| if 
a =1 and viceversa, and s is the index associated with the spin carried by the phonon. 

To perform the numerical calculations one may replace the incoming energies e^, at 
the denominators of Eq. ()34|) with their moduli. The corresponding equation is only 
approximately correct and coincides with the one obtained within the Bloch-Horowitz 
perturbation theory using the Hamitonian in Eq. ()34p |26j . With this replacement 
though, the denominators in Eq. (jMj) never vanish and the framework is therefore well 
suited to obtain numerical results. The two formulas on the other hand necessarily lead 
to very close predictions for since they coincide at the Fermi surface, e^, = 0, where 
pairing is most important. 

In the case of a frequency independent interaction Eq. (jH^ has a simple solution 
also including the level width. This case was treated in Ref . . For instance neglecting 
the phonon-induced interaction in Eq. plj) one obtains 



2e,- 



■/i(e^/,7^/,T), 



where 



n=0 



2e„ 



(35) 



_{hujn + 'y„) + el 

As a first refinement of the theory one can combine Eq. ()35|) and the Bloch-Horowitz 
version of Eq. (jH^ to get 



2e,/ 



duo 



TC 



UJ 



[l-2n^(e,0] 



For a spherically symmetric system the gap, as well as all the single particle properties, 
are independent of the magnetic quantum number m and the equation can be further 
simplified. 



2.5. Multipolar expansion 

Due to the spherical symmetry, the Hartree wavefunctions take the form 0i^(r) = 
Rni{f)Yim{^) ■ Substituting in addition the multipolar expansion for the response 
function (see Eq. (jA.6|l ) the angular integrals in Eq. (j20|l become 



dnYi^{n)Y;^,{n)YLMm 



X 



-1 



' {lm\\YLM\\l'm') 



I I' 



m m' M 

and the same result holds for its complex conjugate. Here 

■(2/ + l)(2L + l)(2/' + l) 



(36) 



(ImWYLMWl'm) 



1/2 



/ /' L 
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is the reduced matrix element and 

/ /' L 

m m' M 

is the Wigner 3-j symbol. 
Recalling that 



(37) 



E 

mm' M 



I V L 
m m' M 



{31 



the gap equation Eq. (jHBj) can finally be written in the rotationally invariant form 



An'l' 121' +1 r (00) . 

n'l' 

1 - 2nF{en'i')]y 



^nLn'l'i^) 



(39) 



'^nZn'«'('^) is matrix element coupled to zero angular momentum and is given by 

-.2 1'+^ 



(00) / \ 
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V2rTTV2i 



L=\l'-l\ 

X Im{ni ^^{r, r', u) - 2Ul nir, r, u)}, (40) 

with fii'L{r) = {lm\\YiM\\l''m')Rn'i'{r)Rni{r). An analogous definition holds for g^i%i,- 

The same multipolar expansion can also be applied to the calculation of the self- 
energy. Eq. (fT^ then becomes: 



21' + 1 



X 



1 n dE (00) 
1 + nB{E) - npiiri'v) , nB{E) + npi^u'i') 



21 + 



+ 



iUn - h ^in'V - E iUn- h ^^n'V + E 



(41) 



where here 



(J 



(00) 



2 

, drr"^ dr' r''^ fii'i(r)fii'i(r') 



X Im{2nL ||(r, r', w) - Ul ti(r, r, cj) - xl ui^, r, uj)}. 

(42) 

Eq. (j39|l is our reference equation for the evaluation of Tc. It is an eigenvalue 
equation for the eigenvector A„;. The critical temperature is the highest temperature 
at which the operator on the RHS has eigenvalue 1. The equation contains a sum 
over n' which needs special care. The contribution of the bare interaction, in fact, 
leads to an ultra-violet divergence, since it arises from a contact approximation to 
the true interatomic potential. Several approaches have been developed to regularize 
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this divergence. They involve a renormahzation of the couphng constant g, which 
is ehminated in favor of the low energy t-matrix {t = ATrh'^ap/m) by means of 
the Lippmann-Schwinger equation, or the introduction of a contact pseudopotential 
g6{r)[{d/dr)r]- |2Z1 I2H]- For a trapped gas, for which only the bare interaction is 
considered, the results of these calculations may be reproduced by directly replacing 
g with the t-matrix, provided one simultaneously introduces a cut-off in the sum at 
^n'l' = [221 I2n]- Because the main scope of the present work is on the phonon- 
induced effects, we treat the bare interaction at this level of approximation. Due to the 
explicit dependence on e„/;/ instead, the phonon-induced interaction has a natural cut- 
off. Numerical calculations have shown that the terms with |e„/i'| ^ give negligible 
contributions. 

3. Results 

In this Section we present the results of our numerical calculations on the quasi-particle 
properties and the critical temperature. As we mentioned in the introduction, the 
calculations were carried out for two specific spherically symmetric infinite 

square well, and a 3D isotropic harmonic oscillator. The infinite square well case is 
interesting for two reasons. First of all, it has very similar properties to those of an 
infinite homogeneous system, and the results of numerical calculations can therefore be 
compared with the predictions of the analytical or semi-analytical formulas. Second, a 
detailed comparison with the harmonically trapped case allows to extract the important 
physical consequences of the discrete shell structure induced by the harmonic trapping. 
The results are shown first for the spherical well, and then for the harmonically confined 
system. In both cases the calculations were carried out with a number of particles fixed 
to = 1000 and three different values of the interaction strength A defined in the 
introduction. 

3.1. Infinite square well and uniform system 

We here consider a system of 1000 particles in a spherically symmetric infinite square 
well of radius R. For these conditions we find that although the BCS coherence length 
(related to the T = pairing gap A(0) by ,^ = h^kp /m'Kls.if))) is generally larger than the 
size of the system, the finite size effects are negligible, and we are thus allowed to relate 
the properties of the system to those of an infinite homogeneous one with the same 
average density. The calculations were carried out for the three different interaction 
strengths A = 0.2, A = 0.4 and A = 0.6, with the temperature fixed at the calculated 
accounting for the direct contact interaction alone. 

Before presenting the results we briefly review the well-known properties of a 
homogeneous system. Ignoring self-energy and induced interaction, the (BCS) gap 
equation leads to the analytical result for the critical temperature 

ksT^o = —-^ep e-'/\ (43) 
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C ~ 0.577 being Euler's constant. 

Also in the weak coupling limit, however, the induced interaction must be included 
in the calculation, since it leads to a constant decrease in the critical temperature by a 
factor (46)^/^ ^2.2 The correct critical temperature for a uniform system in the 
weak-coupling regime is therefore 

„c /o\ 7/3 

This result may be obtained using Eq. ()36|1 by replacing in the expression for the effective 
interaction lip and 11^^ with the unperturbed response function 11° (an approximation 
valid in the weak-coupling limit). Within the present context, one should also ignore 
the self-energy effects {ji,, Re(Sj,) and Z^), and perform an average at the Fermi surface 
of the induced interaction ^T] . 

As the interaction strength A is increased, both the renomalization of the quasi- 
particle properties due to self-energy effects and the collectivity of the modes become 
important. One must therefore keep all terms in Eq. ()36p and solve the full eigenvalue 
equation. For an infinite homogeneous system the problem was considered in Ref. [32] . 
The author used a different set of approximations as compared to the ones we used. He 
considered the renormalization of the quasi-particle properties only at the Fermi surface, 
where he performed an angular average in the spirit of Gorkov and Melik-Barkhudarov 
|30j . At the same time the author kept the full a;-dependence in the quantities, i.e. he 
did not use the one pole approximation. Below we discuss the comparison between our 
results and those of Ref. [32] for the critical temperature. 

Before doing so, we turn to the analysis of the self-energy effects and the quasi- 
particle properties. In Fig. [S] we show the width •y^, of the level z/, calculated using 
Eq. ((SI), as a function of the Hartree quasi-particle energy ^i^. The three curves 
correspond to the values A = 0.2, A = 0.4 and A = 0.6 of the interaction strength. 
The figure shows a very strong dependence of 'y^, on the interaction. In particular there 
is a very large increase in the level width in going from A = 0.4 to A = 0.6. The large 
values of ji, for A = 0.6 question the validity of the one pole approximation for this or 
higher values of the interaction parameter. As we shall see we reach similar conclusions 
also for the harmonically trapped case. 

In Fig. [3 we show the same curve as the one in Fig. [HI for A = 0.4 and compare 
it with that obtained for an infinite system at the same interaction strength using the 
plane waves basis. The calculation for the square well shows a staggering due to the 
discrete shell structure considered. The overall behavior is, nevertheless, very close to 
that of an infinite system. 

In Figs. [SI and [HI we display the real part of the self-energy. We observe that this 
quantity is almost symmetric about the Fermi surface, where it changes sign. The net 
effect is an increase in the density of states at the Fermi level (cf. e. g. [31] and refs. 
therein) . 

We now turn to the results for T^. In Fig. [Tjwe show the value of as a function of 
A. The dashed line corresponds to calculations including the bare interaction alone, the 
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solid line to calculations including all quasi-particle renormalization effects. In Table 
we report the results of the calculations when each single renormalization contribution 
is included. The real part of the self-energy gives rise to an increase in the level density 
around the Fermi surface (see Fig. El) and a consequent increase in Tc due to the 
resulting enhanced effective coupling strength. Both Ziy and •ju decrease the value of 
Tc. The former causes a weaker effective interaction between quasi-particles, as can be 
deduced from Eq. recalling that < 1. Because 7^ measures the energy range 
over which the quasi-particle state is spread due to the coupling to vibrations (lifetime), 
its presence effectively inhibits pairing between quasi-particles [22] • The net effect of 
the induced interaction is to decrease the value of Tc, as can be seen by looking at the 
circles in Fig. |H| There we show the value of for A = 0.2, calculated including only 
the contribution of the effective interaction, as a function of the maximum multipolarity 
included in the calculation of the correlation functions in Eq. pn|l (see the Appendix 
for details on the multipolar expansion). If the overall effect of the combination of spin 
singlet (density) and triplet (spin) modes is repulsive, isolating the two contributions 
shows that density modes provide an attractive effective interaction (see squares in Fig. 
E}, and spin modes a repulsive one (triangles). 

Our results indicate that A = 0.2 does not belong entirely to the weak-coupling 
regime. If this were the case, in fact, inclusion of the induced interaction would result in 
a reduction in Tc by a factor of the order of 2 as predicted in Ref . |HD! ■ This is however 
only true if II'" ~ 11°. When A = 0.2 we indeed find a reduction by factor ~ 2 but only 
when n° is used instead of the full 11^ in the matrix elements in Eq. ()20|1 . Using the 
full we found a reduction only by a factor ~ 1.4. We are lead to conclude that, at 
this density, the collectivity of the modes still plays a significant role. 

For what concerns the comparison with the results of Ref. [32] , we find a reduction 
in Tc in the stronger coupling regimes which is not at all as large as the one reported 
there. We have carefully analyzed our set of approximations and we know that they 
are less and less reliable as the interaction becomes stronger. We are not aware of 
such a careful analysis for the framework of Ref. |H2j which was, on the other hand, a 
pioneering calculation in this field. It is therefore not clear as yet which is the best set 
of approximations to be used, and thus the appropriate values for Tc one should expect, 
as these appear to be strongly model dependent. 

3.2. Trapped system 

Once again the calculations were carried out for the three different interaction strengths 
A = 0.2, A = 0.4 and A = 0.6, with the temperature fixed at the Tc for the system 
interacting via the direct contact interaction alone. 

In Figs. Uni and [TT] we show the width 7,^ and the real-part of the self energy (^,^) 
respectively, as a function of ^j,. Like for a uniform system from Fig. ^Jwe see that 
the levels acquire an increasingly larger width as the interaction strength is increased. 
For A = 0.2 and A = 0.4 the width of the levels at the Fermi surface (^^^ = 0) is small 
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Figure 3. The figure sliows tiie widtli ji, of tiie level v as a function of the hartree 
energy for a spherically symmetric infinite square well. The three curves correspond 
to the interaction strengths A = 0.2 (solid line), A = 0.4 (long dashed line) and A ~ 0.6 
(short dashed line). The energies are measured in the units of /2maR'^ , with R being 
the radius of the well. 
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Figure 4. The curve in Fig. |21 for A = 0.4 (solid line) is compared with the one 
obtained for an infinite homogeneous system with the same interaction strength using 
the plane waves basis (dashed line) . For the comparison the energies are here measured 
in units of the Fermi energy. 



and the one pole approximation is expected to be satisfactory. For A = 0.6 one has 
'-fiy ~ 0.25/ia;o. Since the width is on the one hand much larger than ReEP'^(ej,) and on 
the other almost of the order of the discrete level spacing Huq, the one pole approximation 
is again probably not very meaningful at this or higher values of A, consistently with 
the considerations made for a uniform system. 

The renormalization factor Z^, requires special care when it is evaluated for a finite 
sized system. In Fig. Elwe show ReS^'^(a;) for the level at the Fermi surface. While 
the general trend of the function is definite and follows a curve which resembles that 
of a uniform system, the discrete level structure causes a staggering which makes the 
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Figure 5. The figure shows the real part of the self-energy ReS,y(ej^) of the level v as 
a function of the hartree quasi-particle energy for a spherically symmetric infinite 
square well. The three curves correspond to the same interaction strengths as in Fig. 
Inland the same line codes and energy units are used. 




-0.8 -0.6 -0.4 -0.2 02 0.4 0.6 0.8 

Figure 6. Here we compare the curve in Fig. |S1 for A = 0.4 (solid line) with the 
one calculated for an infinite homogeneous system using the plane waves basis (dashed 
line). The energies are measured in units of the Fermi energy. 





A=0.2 


A=0.4 


A=0.6 




0.86 


0.81 


0.94 




1.05 


1.08 


1.10 


.9+ Z 


0.89 


0.84 


0.80 


9+ 7 


0.97 


0.92 


0.70 


V''ff'+ ReT.P'' + 7 + Z 


0.77 


0.65 


0.49 



Table 1. The table gives the ratio Tc/T^'^™ associated with the indicated contribution 
for the case of a spherically symmetric infinite square well. T^^'° is calculated including 
only g and we find the following results: for A = 0.2, kBT^'""' = 7.3h'^/2maR^, for 
A = 0.4, fcsT]'^" = 19.2fiV2mai?2 ^nd for A = 0.6, keT^''"' = Ah.lh^ /2maR'^ ■ 
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Figure 7. The critical temperature Tc is plotted against the interaction strength for 
the gas in the spherically symmetric inifinite square well. The solid line is the result of 
the calculation made including the full quasi-particle renormalization properties. The 
dashed line corresponds to the calculation using the bare interaction alone. Tc is in 
units of /2maR^ ■ 




Figure 8. Critical temperature of the gas in a spherical symmetric infinite square 
well, with A — 0.2, as calculated from Eq. (|39|l including only the effective interaction 
contribution, as a function of the maximum multipolarity L included in the correlation 
functions appearing in the matrix elements in Eq. H2U|) . L = \ corresponds to the bare 
interaction alone. Convergence is attained around L = 35. Tc is in units of It? /2raaB?' , 
R being the radius of the well. Circles represent the result of the calculation which 
includes the full induced interaction. Squares are obtained accounting only for the 
exchange of density modes and triangles only accounting for spin modes. 



derivative dIieTPj^{uj) / duo not well defined. We have therefore interpolated ReS^^(c(j) 
with a polynomial before taking the derivative. 

For A = 0.2 and A = 0.4 the value of about the Fermi surface has shown a very 
weak dependence on the degree of the polynomial chosen for the interpolation (cf. also 
Ref. jSlj)- It is remarkable that the function ReE^''(ci;) also has a very weak dependence 
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Figure 9. Real part of the phonon self-energy 'ReTjP^{oj) relative to the state ly at the 
Fermi surface (dots), as a function of oj — ep being ep the Fermi energy, at A = 0.2. 
The continuous line is the polynomial curve used to take the derivative and calculate 



on V. The one reproduced in Fig. IHl very closely resembles those evaluated for different 
values of v. Since the staggering increases in size as one moves away from the Fermi 
surface, both the interpolation and the derivative become more and more ill-defined. 
However the contribution of these levels to the gap equation is negligible, for this reason 
the procedure used has demonstrated to be solid and so has the calculated value of T^. 
For A = 0.6 on the other hand the staggering is large also at the Fermi surface. This 
once again points at the fact that the approximation used may not be entirely applicable 
for this interaction strength. 

We now turn to discuss the effects on Tc. In Fig. El we show the results for 
Tc as a function of A. The two curves correspond to Tf. as evaluated using the bare 
interaction alone (dashed curve), and including all quasi-particle renormalization effects 
(solid curve). 

General features - In Table |2l we report the results of the calculations when the 
specific effects of each renormalized quasi-particle property is isolated. The effects of 
the self-energy renormalizations are analogous to those we dealt with in the uniform 
case. The real part ReS^'^ causes an increase in the critical temperature since, as 
before, it leads to an effective increase of the density of levels at the Fermi energy. On 
the other hand, both Zy and depress the value of Tc, the former by causing a weaker 
effective interaction between quasi-particles and the latter by inhibiting pairing between 
quasi-particles [29] . 

Specific features - In the weak coupling regime A = 0.2 the critical temperature is 
hardly affected by the renormalization of the quasi-particle properties. While this is 
expected as far as 71,, Zy and ReS^'^ are concerned, it is instead surprising as regards 
the effect of the induced interaction. For a uniform system at A = 0.2 we have seen 
that, if the induced interaction is calculated using only 11° one finds a decrease in Tc by 
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A=0.2 


A=0.4 


A=0.6 




0.92 


0.86 


0.90 




1.01 


1.03 


1.03 


9+ Z 


0.98 


0.90 


0.93 


5+ 7 


0.97 


0.85 


0.76 


V''ff+ ReY.P^ + 7 + Z 


0.88 


0.66 


0.61 



Table 2. The table gives the ratio T^/T^'^^'^ associated with the indicated contribution 
for the case of a sphericaUy symmetric harmonic potential. T^'*"'® is calculated including 
only g and we find the following results: for A 0.2, fcsT^'*™ = 5.03 x 10~^fia'o, for 
A = 0.4, /cbT^^"° = 0.34fta;o and for A = 0.6, fcsT^^^'^ = 0.75ftu;o- 

a factor of the order of 2. This effect is reduced, but still present, when the collectivity 
of the modes is accounted for. Almost no reduction is obtained instead for a trapped 
gas, whether one uses the lowest-order correlation function 11° or one accounts for the 
full collectivity of the modes. The absence of the reduction descends here from the 
discrete level structure of the harmonic trap. This can be understood by comparing the 
strength functions of a uniform system with those of a harmonically trapped one for 
the same interaction strength, see Figs. lAlf IA2| IA5I and IA6I The largest contributions 
to the sum in Eq. ()36|) comes from the condition = e^i = uj = 0, which corresponds 
to particles at the Fermi surface and to the exchange of a phonon with zero frequency. 
The discrete level structure causes the values of the strength functions to be negligible 
near u; = 0, and consequently the effect of the phono n-induced interaction to be small. 
As the strength of the interaction increases (A = 0.4 and A = 0.6) the renormalization 
of the quasi-particle properties has a larger and larger effect on T^. All the terms are 
increasingly more important but it is the quasi-particle width that dominates over the 
others, as can be seen in Table El We recall however that the case A = 0.6 is beyond 
the validity of the one pole approximation. 

4. Conclusions 

In this paper we have calculated the renormalization of the quasi-particle properties 
in an interacting Fermi gas due to the emission and reabsorption of density and 
spin fluctuations of the system. In particular we have looked, within the one pole 
approximation, at self-energy (S) and screened interaction {V^^^) effects. The real part 
of S leads to a shift in the energy of the quasi-particles and to a reduction in the 
strength of the quasi-particle peak of the associated propagator (by a factor Z < 1). At 
the same time a non-vanishing imaginary part of S results in a non-vanishing width 7 
of the levels, or equivalently in a finite lifetime of the quasi-particles. The exchange of 
collective modes between quasi-particles leads to an effective interaction which has an 
overall repulsive character. This is due to the dominance of the spin triplet (repulsive) 
modes over spin singlet (attractive) modes, and leads to a screening of the bare attractive 
interaction. 
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Figure 10. The width 7^ of the level v is plotted as a function of the hartree energy 
for a spherically symmetric harmonic trap. The three curves correspond to the 
interaction strengths A = 0.2 (solid line), A = 0.4 (dashed hne) and A = 0.6 (dot- 
dashed line). The energies are measured in units of Huiq, with being the oscillator 
frequency. 
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Figure 11. The figure shows the real part of the self-energy ReE^''(ei,) of the level v as 
a function of the hartree quasi-particle energy for a spherically symmetric harmonic 
trap. The three curves correspond to the same interaction strengths as in Fig. 1101 and 
the same line codes and energy units are used. 



Finally we have calculated the effect of these properties on the critical temperature 
to the superfluid state. The calculations were carried out numerically for the two 
illustrative cases of an isotropic infinite square well and for a harmonic trap. In both 
cases the calculations were carried out for a sytem of 1000 particles. We have included 
selectively one or the other of the renormalized quasi-particle properties and found that 
these have different consequences of T^. When all of them are accounted for, the overall 
effect is a reduction of Tc which depends both on the strength of the interaction and 
on the geometry of the system. The results are reported in Sec. El As the strength 
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Figure 12. Critical temperature of a gas of 1000 particles in a spherically symmetric 
harmonic oscillator. Tc is in units of huo- The dashed curve corresponds to Tc 
calculated using the bare interaction alone, the solid one is obtained including all 
quasi-particle renormalization effects. 
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Figure 13. Critical temperature of the harmonically trapped gas with A = 0.2, 
as calculated from Eq. H39|) including only the effective interaction contribution, 
as a function of the maximum multipolarity L included in the correlation functions 
appearing in the matrix elements in Eq. (|20|l . L = 1 corresponds to the bare interaction 
alone. Convergence is reached around L = 27. Tc is in units huio, with ojq being the 
trap frequency. Circles represent the result of the calculation which includes the full 
induced interaction. Squares are obtained accounting only for the exchange of density 
modes and triangles only accounting for spin modes. 

of the interaction increases, although the one pole approximation becomes less and less 
reliable, our calculations clearly indicate that the effect of fluctuations becomes more 
and more important and modifies the critical temperature, as compared to the one 
obtained using the bare quasi-particle properties, in an important way. This result is 
very relevant for current experiments with atomic Fermi gases. In particular it shows 
that the standard theories of the crossover from BCS to BEC, which are currently used to 
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predict the value of the critical temperature also in the strongly interacting regime may 
be wrong by a significant amount as the interaction becomes stronger. Our approach 
is well suited to calculate the first corrections to the quasi-particle properties and to 
the critical temperature which arise when one moves towards the strongly interacting 
regime originating from the BCS weak coupling side. 



Appendix 

Appendix A.l. The correlation functions 
Define: 

Yi^p{x,x') = -e^^Tr{pr,[n„(x)n^(x')]} 

where h^{x) = h^{x) — {'ijjl{r)ijj„{r)), with h^{x) = ■ijj1{x)'ijj(j{x) and x = (r, r). The 
indices a, jS and a can assume the two values f and |, the operator p is given by 
p — exp[H^ — /xiV], with N — J d^r '?/'J.(r)^o^(r) and Tt is the ordering operator 
relative to the imaginary time r. Finally exp(— is the grand partition function of 
the system, and Q the thermodynamic potential. 

Density correlation function - Let us now introduce the density correlation function 
with the definition: 

Xlp{x, x') = -e^"Tr{p Tr[n{x)h{x')]} (A.l) 

where n — + h]^. This correlation function is related to the previously defined as 
follows: 

n,(x,x') = 2nT^(x,x') + 2n^i(a;,a;'), (A.2) 

where we have used n||(a;, a;') = n^(a;, a;') and Ii^^{x,x') = Ilii{x,x'). 

Spin correlation function - The spin correlation function has components along the 
three directions x, y and z, defined as follows: 

U„X^,x') - e^''Tr{pT,[ai{x)ai{x')]} (A.3) 

where (jj = "^^f^ i^ti^) i^i) a^i^pi^) ^ and {<Ti)ap being the Pauli spin matrix element. The 
spin correlation function in the z direction is related to IIq,^ by: 

U,^{x,x') = 2n^^{x,x')-2U^^{x,x'). (A.4) 

Finally since the choice of the orientation of the axes is arbitrary we must also have 
Ila^{x,x') = Il(jy{x,x') — Ila^{x,x'). The Ua/s therefore contains all the information 
needed. 

Ha/B admits the Lehmann representation: 
U^,{vy,iu;.)^e^%-^-^-e-^--) ^ ^"^-'""^^^'^f f 

,„ [-t^n -t^m) 

where the complete basis of eigenstates \^n) introduced above is in principle the exact 
basis of excited states of the system. Below we shall assume specific approximations for 
the |$n)'s. 
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Appendix A. 2. The RPA approximation 

To evaluate one has to introduce specific approximations. The lowest order 
approximation is the Hartree-Fock particle-hole one, in which case: 

^ V Xo(r,r',?CL;„) j 

As is well known a better approximation for is the RPA self-consistent one. In this 
other case Hap is assumed to satisfy the equation (for each given ia;„): 

nf^(r,r') nf^(r,rO \ ( Xo(r,r') 
nj^(r,r') nj^(r,r') ) \ xo(r,r') 

-|- J d'^rid'^r2 

g5(r^-r2) \ ( n^^(r2,r') n^^(r2,rO 
g5{r^,V2) j \ n\^{v2,r') nj^(r2,r') 

In compact notation one can write W^.^ — 11^ = Xo + 9^Xo + ■■■ ^iid n^j^ = nj| = 
yXo + 9\t + - Consequently = 2xo/[l - «/Xo] and W^^ = 2xo/[l + ^Xo]- 

Appendix A. 3. Multipolar expansion 

The spherical symmetry of our system allows expansion in multipoles: 

n„^(r,r') = Y,[^c.p{ry)]LYLM{miM{^') (A.6) 

LM 

where Q, and Q! are the angular coordinates of r and r' respectively. 
Appendix A. 4- Strength function and sum rules 

An essential tool to study the collective modes is the strength function. This is defined 
as 

S{F, E) = e^"" ^(e-'^^" - e'^"^-) 

nm 

X \{m\F\n)\^5{E-Er,r^), (A.7) 

where F = J d^rF^{r)ipl{r)tl)^{r) is the operator which corresponds to external fields 
F|(r) acting on the t particles and -F|,(r) acting on the J, ones. Enm stands for En — Em- 
The strength function is related to the retarded density correlation function by 

S{F,E) = ~Y. j dhd'T'F;{v)F„,{v')lm{Iil^,{vy,E)}. 

Density modes (i.e. with spin s = 0) are excited by external fields such that 
F|-(r) = -^^(1"), while spin modes (with s — 1,8^ — 0) are excited by fields with 
F^(r) = -F^(r). 
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Figure Al. The figure sliows tlie strengtli functions for density and spin modes with 
multipolarity L = 2, for a system of 1000 particles in a sperically symmetric harmonic 
trap and interaction strength A = 0.2. The solid line shows the strength functions 
calculated with the Hartree-Fock correlation function II''^, in this case density and 
spin response coincide. In the RPA the two responses and distinct and are given by 
the dashed line and dot-dashed line respectively. The energy on the x-axis is measured 
in units of the oscillator level spacing Hluq. 
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Figure A2. Same as Fig. lAll for the L = 12 mode. 



In the particular case of a density mode of given multipolarity L, for which 
F|(r) = -F|(r) = r^Yijifi^), upon use of Eq. ()A.6|) . this becomes 

Sl{E) = -^Y1 j drdr'r'^+\'''^hm{[UlAr,r,E)]L}. 

era' 

Also useful is the energy weighted sum-rule jHHj 

J dES{F,E)E = ^Y. j d''r\VF^{v)\^n^{r) (A.8) 
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Figure A3. Same as Fig. lAll for the interaction A = 0.4. 




Figure A4. Same as Fig. lAll for the L = 12 mode and interaction A — 0.4. 

When F{r) = t^YlmI^) it takes the simple form 

J dE Sl{E)E = -^J2j dr' r^+h'^^^ x Im{[nO^,(r, r', E)]^} 

(TO-' 

= L(2L + 1)/ drr^^nJr). 

2m Jo 

and can be used as a check for the numerical calculations. 



(A.9) 
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